30, 2010 0:23 WSPC/INSTRUCTION FILE DTK 



International Journal of Modern Physics D 
© World Scientific Publishing Company 



Coherent versus coincidence detection of gravitational wave signals 
from compact inspiraling binaries 

S. Dhurandhar 1 , H. Mukhopadhyay 1 , H. Tagoshi 2 and N. Kanda 3 

1 IUCAA, Postbag 4, Ganeshkind, Pune -411 007, India. 
2 Department of Earth and Space Science, Graduate School of Science, Osaka University, Osaka 

560-0043, Japan. 

3 Department of Physics, Graduate School of Science, Osaka City University, Osaka 558-8585, 

Japan. 

We compare two multi-detector detection strategies, namely, the coincidence and 
the coherent, for the detection of spinless inspiraling compact binary gravitational wave 
(GW) signals. The coincident strategy treats the detectors as if they are isolated - com- 
pares individual detector statistics with their respective thresholds while the coherent 
strategy combines the detector network data phase coherently to obtain a single detection 
statistic which is then compared with a single threshold. In the case of geographically 
separated detectors, we also consider an enhanced coincidence strategy because the usual 
(naive) coincidence strategy yields poor results for misaligned detectors. For simplicity, 
we consider detector pairs having the same power spectral density of noise, as that of 
initial LIGO and also assume the noise to be stationary and Gaussian. We compare the 
performances of the methods by plotting the receiver operating characteristic (ROC) for 
the two strategies. A single astrophysical source as well as a distribution of sources is 
considered. We find that the coherent strategy performs better than the two coincident 
strategies under the assumptions of stationary Gaussian detector noise. 

1. Introduction 

Inspiraling binaries are one of the most promising sources for the first detection of 
gravitational waves (GW). The post Newtonian approximation methods accurately 
describe the phasing of the waveform - about a cycle in a wave train ~ 10 4 cycles 
long. This makes it amenable for matched filtering analysis. The best available esti- 
mates suggest that the expected number of neutron star (NS)-NS binary coalescence 
seen per year by ground based interferometers is 7.1 x 10 -3 — 0.12 for initial detec- 
tors and 38 — 6.6 x 10 2 for advanced detectors^. In recent years, a number of ground 
based detectors are producing sufficiently interesting sensitive data and analysis of 
network data is highly advisable. The advantages of multi-detector search for the 
binary inspiral is that, not only does it improve the confidence of detection, it also 
provides directional and polarisational information about the GW source. 

Two strategies currently exist in searching for inspiraling binary sources with a 
network of detectors: the coherent and the coincident. The coherent strategy involves 
combining data from different detectors phase coherently, appropriately correcting 
for time-delays and polarization phases and obtaining a single statistic for the full 
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network, that is optimized in the maximum likelihood sense. On the other hand, 
the coincident strategy matches the candidate event lists of individual detectors 
for consistency of the estimated parameters of the GW signal. However, the phase 
information is ignored and as also the detectors are considered in isolation. 

The question arises as to which strategy performs be tter . On simple waveforms 
the analysis has been performed by Finn and Arnaud et aP^l Both these works have 
shown that the coherent strategy performs better than the coincident strategy. We 
consider here the astrophysically important source, namely, the inspiraling binary 
and report main ly the results of our work which has been described in detail in 
the papers I 4 | 5 | 6 | CO mpare the strategies by plotting the Receiver Operating 
Characteristic (ROC) curves, which is the plot of detection efficiency versus the 
false alarm rate. We broadly consider the two cases of co-located aligned detectors 
and geographically separated misaligned detectors. In the co-located case we further 
consider two subcases of (i) uncorrelated noise, and (ii) correlated noise. 

2. The correlation statistic 

For the inspiraling binary, in the Fourier domain we assume the spinless restricted 
post-Newtonian (PN) waveform h 1 (/) at detector /: 

h\f) = ME'r 7 ' 6 exp »[¥(/; t c , 5 C , r , r 3 ) + 2 7 r/At / ] , (1) 

where we take the phase \& given by the 3 PN formula. The extended beam pat- 
tern functions ^ E 1 encode the orientation and direction parameters of the source 
and detectors; the parameters t c ,S c , are respectively, the time of coalescence, phase 
of coalescence; the quantities tq , T3 are the chirp time parameters which are inde- 
pendent functions of the two masses of the stars comprising the binary, the At 1 
denote time-delays at the detector / with respect to a fiducial detector and J\f is 
the amplitude depending on the masses and the distance to the source. 

In matched filtering, it is natural to define the scalar product (a, b) of two real 
functions by, 

W) = ! f>w™, (2) 

where, we use the Hermitian property of Fourier transforms of real functions. Sh(f) 
is the one sided power spectral density (PSD) of the noise which is assumed to be 
the same for all the detectors. Then the normalised templates sq, s t / 2 corresponding 
to the intrinsic parameters p = {tq, T3} are defined via the equation, 

h T (f; pi, t c , S c ) = A I (s (f; Pi, tc) cos 5 C + s 7T / 2 (f; Pi, t c ) sin S c ) . (3) 

We have the relation s n / 2 (f] Pi, t c ) = i§o{f; pi, t c ) where jU, represents a grid point 
in the intrinsic parameter space and the normalisation of templates requires that the 
scalar products (s , Sq) — {s^n, s^/2) = 1. We then define the complex correlation 
C 1 = Cq + ic^/2 where, the real correlations Cq and c^./ 2 are obtained by taking 
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the scalar products of the data x 1 with sq and s^/ 2 respectively; that is, Cq ff , 2 = 
( s o,ir/2, xI )- Then the single detector statistic for detector I is just A/ = \C I \ 2 . To 
decide detection in a given single detector I, A/ is maximised over the template 
parameters pii,t c and compared with a preassigned threshold. 

3. Coincidence versus coherent detection 

For coincidence detection, the procedure is as follows: 

• Choose the same threshold A* for the two detectors. 

• Prepare two candidate event lists such that A/ > A*, I = 1,2. Look for 
pairs of candidate events, each candidate event coming from a different 
list, such that the sets of estimated parameters £ c ,to,T3 match - |Ai c | < 
Ayyic, |At | < AyyT , |Ar 3 | < AyvT3 where the AA denotes the difference 
in the measured parameter A, where A stands for any of the parameters 
tcTQ^T^. The allowed error box is denoted by AyyA. A quadratic sum of 
the noises is taken to determine the error box. We fix this box by performing 
simulations, so that the final probability of not losing an event is 0.97 - on 
each parameter, we atmost allow 1 % loss in events. For geographically 
separated detectors, for fixing the window size in t c , the light travel time 
between the detectors is taken into account and added in quadratures to 
the errors due to noise in each detector. 

On the other hand, coherent detection involves combining data streams in a 
phase coherent manner so as to effectively construct a single, more sensitive detector. 
For the case of two misaligned detectors the network statistic is given by^l 

A = ||C|| 2 = IC 1 ! 2 + |C 2 | 2 = {elf + (cl /2 f + (c 2 ) 2 + (c 2 /2 ) 2 , (4) 

where C 1 is the complex correlation of the 7-th detector (1=1,2). For two aligned 
colocated detectors the statistic is different ^ and is given by: 

A = i|C 1 + C 2 | 2 . (5) 

For aligned colocated detectors we consider two subcases: (i) uncorrelated noise, (ii) 
correlated noise. 

4. Results 

4.1. Colocated aligned detectors 

We first consider the case of co-located aligned detectors. This case would be of 
significance to the two LIGO detectors at Hanford and other similar topologies 
envisaged elsewhere in the future such as LCGT. We compare the performances 
of the two strategies by plotting the ROC curves for uncorrelated noise (left) and 
correlated noise (right) in Fig. (JTJ). For plotting these curves the false alarm and 
the detection probabilities must be computed as a function of the threshold A* and 
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then plotted versus each other parametrically by varying the parameter A*. The 
correlation parameter eo is taken as the weighted average of a frequency dependent 
correlation e(/) where (ni(/)n|(/')) = ^e(f)S h (f)S{f - /'), ni(t) being the noise 
in the I th detector. Details may be found inSEl. While performing the simulations 
there are many subtleties such as estimating the number of independent templates, 
the error window size computation etc. the discussion of which we have omitted 
here although it is nonetheless important. 




Fig. 1. The ROC curves for single detector, coincident and coherent analysis for an injected 
signal of SNR=10. The solid line, dashed line and the dotted line correspond to the theoretical 
ROC curves for coherent, coincident and single detectors respectively. The left hand side figure 
corresponds to uncorrelated noise while the right hand one to correlated noise with eo = 0.3. 

It is clear from the figures that the coherent strategy is far superior in this case. 
4.2. Geographically separated misaligned detectors 

We now consider the case of geographically separated detectors which are then also 
normally misaligned. When the detectors are misaligned, the sky coverage for the 
usual coincidence detection is very poor which leads to intolerable false dismissal. 
Thus another coincidence strategy is devised which we call enhanced coincidence. 
The usual coincidence strategy we then call naive coincidence. Enhanced coincidence 
strategy is formulated as follows: 

• Choose a low threshold Aq (we choose Aq ~ 16), and prepare two candidate 
event lists such that A/ > AJJ, I = 1,2. 

• Look for a pair of candidate events, the events coming from separate lists, 
such that the sets of estimated parameters match within the error- window. 
The procedure is the same as the co-located case, except for the parameter 
t c in which the distance between the detectors enters. 

• Choose the final (high) threshold A* > 2Aq and construct the final statistic 
A = Ai + A 2 and register detection if A > A*. 

Note that although this statistic looks formally like the coherent case, the mass 
parameters for the templates in the two detectors do not have to be the same 
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Fig. 2. The ROC curves have been plotted for coherent, naive coincidence and enhanced coinci- 
dence for uniformly distributed sources. The left hand figure is drawn for the LIGO, Livingston 
and Virgo pair of detectors, while the right hand figure corresponds to the LIGO Livingston and 
LIGO Hanford detector pair. 



(they however must be close enough so that they lie in a error window). Thus 
this is not matched filtering while coherent detection is. In this strategy the sky 
coverage is better than naive coincidence. As seen from the Fig. @, coincidence 
strategy performs far better than the naive coincidence strategy, but we see that 
the detection probability for the coherent strategy is still superior by around 5% 
for the same false alarm rate. 



5. Concluding remarks 

Although the coherent strategy is superior to coincident strategies, the difference 
between the coherent and enhanced coincident strategies is small. Only a relative 
improvement of about 5% in the detection probability is obtained with the coherent 
strategy. One may ask whether there is any practical advantage in using the coherent 
strategy in the case of two misaligned detectors. Note however that the coherent 
method is not so computationally expensive compared with two coincident methods, 
since we do not take cross correlation of two detectors' data in the coherent strategy. 
Thus overall, we conclude that the coherent strategy is a good detection method. 

However, the above results assume stationary Gaussian noise. But we know that 
the current real data are neither stationary nor Gaussian. In coincidence detection, 
the requirement of the consistency of estimated parameters in an error window 
acts as a powerful veto ^ to veto out fake events generated from non-Gaussian 
noise. On the ther hand in coherent detection as yet no such obvious veto has 
been developed. This however does not rule out the possibility that a powerful veto 
cannot be constructed for the coherent strategy. In future we propose to work on 
this aspect of the problem. Perhaps a judicious combination of the two methods 
might be an effective way of dealing with this problem. 
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